Continuum-type stability balloon in oscillated granular layers 
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The stability of convection rolls in a fluid heated from below is limited by secondary instabilities, 
including the skew-varicose and crossroll instabilities. We observe a stability boundary defined by 
the same instabilities in stripe patterns in a vertically oscillated granular layer. Molecular dynamics 
simulations show that the mechanism of the skew- varicose instability in granular patterns is similar 
to that in convection. These results suggest that pattern formation in granular media can be 
described by continuum models analogous to those used in fluid systems. 
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A system comprised of a large number of discrete 
grains behaves like a fluid under some conditions and like 
a solid under others, but can also display behavior unique 
to the granular state JjJ . Despite substantial recent inter- 
est in both the statics and dynamics of granular systems, 
there is no unifying theoretical description of granular 
materials, and it has been argued that a general local, 
continuum description of granular media, analogous to 
hydrodynamics, is unlikely to exist 0-0]. In this Letter 
we report on a study of the stability of stripe patterns 
that form when a horizontal layer of granular material is 
oscillated vertically |4|-|(| ■ In both experiments and sim- 
ulations, we observe behavior strikingly similar to that 
seen in fluid dynamical systems j7j. Our results sug- 
gest that a continuum description of pattern formation 
in granular media is possible. 

Recent experiments on pattern formation in oscillated 
granular layers [[|-^| have demonstrated the existence of 
stripe, square and hexagonal patterns pl^J, as well as 
localized structures called oscillons ||, as the frequency 
and amplitude of the vibration and the layer thickness 
are varied. At onset, the patterns are subharmonic, that 
is, they oscillate with a frequency equal to half the driv- 
ing frequency. The patterns are similar in spatial struc- 
ture to those observed in fluid dynamical systems, most 
notably parallel convection rolls in a thin layer of fluid 
heated from below (Rayleigh-Benard convection) [QJ|] 
and standing surface waves in a vertically oscillated liquid 
layer (the Faraday instability) . Molecular dynamics 
simulations Jic|jll|l have reproduced the various granu- 
lar patterns at values of the control parameters equal to 
those used in the experiments JO]. Some aspects of the 
granular patterns have been described using various mod- 
els for the dynamics 0|, amplitude equations flT3| , and 
iterative maps p4| . However, a rigorously derived theo- 
retical description of pattern formation in oscillated gran- 
ular media does not yet exist, and the extent to which 
the well-developed understanding of patterns in fluid sys- 
tems Q can be applied to granular material has not been 
established. 



Our experimental apparatus is similar to that used in 
Refs. @,|. Bronze spheres 150-180 /j,m in diameter are 
contained in a cylindrical cell 14.7 cm in diameter. The 
cell sidewall and top lid are made of Plcxiglas. The alu- 
minum bottom plate is flat out to a radius of 5.0 cm, then 
slopes upwards to the sidewall at an angle of 2°. As do 
fluid convection rolls, the stripes prefer to form perpen- 
dicular to the cell wall; the "beach" reduces the effect of 
the sidewall on their orientation. The cell is evacuated, 
and an electromagnetic shaker, driven sinusoidally, oscil- 
lates the cell vertically at a frequency /. T, the amplitude 
of the acceleration relative to the gravitational accelera- 
tion g, is determined using an accelerometer mounted on 
the bottom of the cell. The pattern is illuminated with 
a ring of LEDs encircling the cell and strobed at f/2. 
Images of the patterns are recorded using a digital video 
camera mounted above the cell. 

Experiments were performed on layers of thickness 
N (scaled by the mean particle diameter) in the range 
6.6 < N < 11.3 , in a frequency range over which stripe 
patterns formed. The layer was flat for V < T c w 2.5, 
at which point stripes appeared via a subcritical bifur- 
cation. Typically, T was increased suddenly from below 
onset to a value in the stripe regime. We observed time- 
dependent patterns which were predominantly stripe- 
like but with curvature, and which contained point de- 
fects, grain boundaries, spiral defects Q, and sidewall 
foci. Stripes persisted as T was increased, until the sys- 
tem (for most conditions) underwent a subcritical tran- 
sition to hexagons 0] . Below this transition we observed 
two secondary instabilities of the stripes — the skew- 
varicose and crossroll instabilities |?],|| — which limited 
the wavenumber range over which stripes were stable at 
a given frequency. Similar phenomena occured in cells 
without a beach. 

Figure [l] shows the localized skew-varicose instability 
observed when the local wave number of the pattern be- 
comes too large. An initially straight pattern of stripes 
develops a distortion. Eventually one or more stripes 
pinch off and dislocation defects form, which then prop- 
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agate away by a combination of climbing and gliding. 
Since one or more stripes are destroyed by this process, it 
leads to a decrease in the local wave number. Analogous 
behavior has been observed in Rayleigh-Benard convec- 
tion P^|l6[ and the Faraday instability |0 . 




FIG. 1. A sequence of images showing the skew-varicose 
instability of a stripe pattern and the formation of a pair 
of defects, (a) Time t — —4.0 s, straight stripes with wave 
number (k) = 10.8 cm -1 ; (b) t — s, skew- varicose distortion 
is visible; (c) t = 0.5 s, defects form; (d) t = 1.0 s, defects 
move away, leaving a stripe pattern with (k) = 9.9 cm - . 
F = 3.10, / = 35 Hz, N = 9.2, and a 4.7 cm square region of 
the cell is shown. 

Figure ^| shows an example of a crossroll instability, 
which is again analogous to that seen in fluid convection 
0,0. Crossrolls are observed to occur via two pro- 
cesses, and, for most conditions, occur when the local 
wave number becomes too small. Squares and/or per- 
pendicular stripes can develop locally, or, as in Fig. |, a 
region of stripes can be invaded by a region of perpen- 
dicularly oriented stripes. 

To determine the wave number k at which the stripe 
pattern becomes unstable to these instabilities, we used 
the method of Ref. [[l9| to determine k at each point 
in a small region of the pattern in which the instability 
occurs, k was averaged over the region of interest for 
a time series of images bracketing the instability. The 
mean wave number (k) is plotted in Fig. [| for the skew- 
varicose event of Fig. |[ (k) increases up to the onset 
of the instability, then suddenly changes when a stripe 
is pinched off and defects form. The standard deviation 
a of (k) increases as the pattern gradually becomes dis- 
torted and is relatively large immediately following the 
instability. After the defects have traveled away, a stripe 
pattern with a smaller wave number remains and a is 
again low. In Fourier-space the skew-varicose instability 
is signalled by the local appearance of Fourier power at 
an oblique angle — typically 15° to 25° — to the wave 
vector of the original stripe pattern. 




FIG. 2. A sequence of images showing the crossroll insta- 
bility of a stripe pattern, (a) Time t — —2.4 s, straight stripes; 
(b) t = s, crossrolls invade from the top edge of the image; 
in the center of the image (k) — 6.8 cm -1 ; (c) t — 3.5 s; (d) 
t — 5.0 s, the pattern is predominantly perpendicular stripes 
with (k) = 7.0 cm" 1 . T = 2.82, / = 26.8 Hz, N = 11.3, and 
the region shown is 4.7 cm square. 
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FIG. 3. The average wave number {k) (top) and standard 
deviation a (bottom) for the skew-varicose event illustrated 
in Fig. h . The times corresponding to the four images shown 
in Fig. I are indicated above the graph. Between the dashed 
lines the pattern is not locally stripe-like due to the presence 
of defects, and the values of (k) and a are plotted as small 
symbols only to illustrate the change that occurs when the 
defects form. 

FigureH shows analogous data for the crossroll event 
of Fig. pT In this case (k) decreases prior to the onset 
of the instability. When the crossrolls invade, the lo- 
cal pattern is disrupted and a increases. At the end of 
the time interval shown, the region contains a pattern of 
straight stripes perpendicular to and with k larger than 
the original stripes. In Fourier space, the development of 
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crossrolls is accompanied by the appearance of power at 
90° to the original wave vector. 
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FIG. 4. The average wave number (k) (top) and standard 
deviation a (bottom) for the crossroll event of Fig. ^| The 
times corresponding to the images shown in Fig. |^ are in- 
dicated above the graph. Values of (fc) and a between the 
dashed lines are plotted as small symbols to illustrate the 
qualitative change that occurs when the crossrolls invade. 

By analyzing a number of similar events over a range 
of r, /, and N, we can construct a stability diagram 
for the stripe pattern analogous to the Busse balloon ||] 
for Rayleigh-Benard convection. We take the maximum 
value of (k) before a stripe pinches off as the wave number 
corresponding to the skew-varicose instability, and the 
minimum value of (k) before the appearance of crossrolls 
as that for the crossroll instability. 

Figure || is the resulting stability diagram for N = 9.2 
and / = 29.7 Hz. It indicates that there is a limited range 
of k over which the stripe pattern is stable. The stability 
range is bounded at high k by the skew-varicose insta- 
bility and at low k by the crossroll instability. Similar 
stability boundaries, but with slightly different shapes, 
were measured at a number of different values of N and 
/. However, for the depth in Fig. [s], the boundary is 
different at high frequencies: At / = 40 Hz, the stripes 
are unstable to crossrolls on both the high and low k 
sides of the stability boundary, and at the same depth 
and / = 45 Hz, squares rather than hexagons limit the 
range of existence of stripes at high V. These stability 
boundaries are very similar to those for Rayleigh-Benard 
convection rolls: For low Prandtl number fluids |2Q| ; the 
skew- varicose instability forms the high-fc limit of the sta- 
bility balloon, while crossrolls are observed on the low-fc 
side |^, [l5| , [l6|] . For higher Prandtl number, the crossroll 
instability occurs at both high and low k [| 18 1. 
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FIG. 5. The stability boundary for stripe patterns at 
N = 9.2, / = 29.7 Hz. Squares indicate the crossroll in- 
stability and circles the skew-varicose instability. Each point 
is the average over several events; the error bars are stan- 
dard deviations. The lower dashed line indicates the point at 
which the pattern fills the cell as F is increased, and the upper 
dashed line marks the transition from stripes to hexagons. 

We have also observed skew-varicose and crossroll in- 
stabilities in event-driven molecular dynamics simula- 
tions which have been previously validated against exper- 
iment jnj. We simulated the motions of 60000 particles 
in a square cell of side 100 particle diameters, with pe- 
riodic boundary conditions in both horizontal directions. 
A pattern of straight stripes with k either too large or 
too small to be stable at the frequency of the simula- 
tion was used as the initial condition. Time sequences 
of the resulting instabilities are shown in Fig. ||. As in 
the experiments, the skew- varicose instability leads to a 
decrease in k while the crossroll instability leads to an 
increase in k. 

Patterns of rolls in fluid convection are well described 
by partial differential equations for the amplitude (near 
onset) and the phase of the pattern |/J. The cross- 
roll instability exists in an amplitude equation descrip- 
tion that allows for the existence of two perpendicular 
modes. The skew-varicose instability, however, depends 
on the existence of a large scale vertical vorticity (or mean 
flow) , which couples to roll curvature |2l| . To determine 
whether an analogous vertical vorticity exists for skew- 
varicose instabilities in the granular layer, the particle ve- 
locities from the simulations were averaged spatially over 
depth and temporally over one oscillation of the pattern. 
The resulting two-dimensional velocity field v was then 
low-pass filtered to remove spatial frequencies due to the 
pattern itself. Finally, we calculated the squared vortic- 
ity, ui 2 = (V x v) 2 , and averaged it over the entire cell. 
The result is shown in Fig. ||(c). The vertical vorticity 
increases as the skew- varicose instability develops, then 
sharply decreases after a stripe pinches off. In contrast, 
there is no change in (w 2 ) for the crossroll instability. 
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FIG. 6. Instabilities of a stripe pattern from simulations of 
the oscillating layer with T = 3.0, N = 5.4. (a) Skew-varicose 
instability at f / y g/H = 0.37; kH — 1.7 initially and 
and 1.4 after the instability. (b) Crossroll instability at 
f/y/g/H = 0.40; kH = 1.36 before and 1.7 after. Here 
H is the layer depth and T is the oscillation period, (c) 
Mean-square vertical vorticity as a function of time for the 
instabilities shown in (a) and (b). 

We have shown that the range of stability of stripe 
patterns which form in an oscillated granular layer is lim- 
ited by the skew- varicose and crossroll instabilities, and 
that the stability boundary of stripe patterns in gran- 
ular layers is similar to that of straight rolls in a con- 
verting fluid. Our simulations produce the same insta- 
bilities, and indicate that vertical vorticity plays a role 
in the skew-varicose instability, as it does in fluid sys- 
tems. These results demonstrate a clear correspondence 
between patterns in oscillated granular media and hydro- 
dynamic systems, and suggest that, despite the unique 
physical properties of the oscillated granular layer, a con- 
tinuum description of this system should be possible. 
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